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Abstract 

We consider the nonequilibrium dynamics of a a real scalar field in a degenerate 
double-well potential. The system is prepared in the lowest free energy state in one of 
the wells and the dynamics is driven by the coupling of the field to a thermal bath. 
Using a simple analytical model, based on the subcritical bubbles method, we compute 
the fraction of the total volume which fluctuates to the opposite phase as a function 
of the parameters of the potential. Furthermore, we show how complete phase mixing, 
i.e. symmetry restoration, is related to percolation, which is dynamically driven by 
domain instability. Our method describes quantitatively recent results obtained by 
numerical simulations, and is applicable to systems in the Ising universality class. 
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It is well-known that an Ising ferromagnet loses its net magnetization above a certain 
temperature. The simplest mean- field description of this phenomenon is based on the 
Ginzburg-Landau model, in which the average magnetization M = {V^ 1 ) J dV M(x) is 
the order parameter, and the thermodynamic potential is expanded to fourth power in M, 
V(M) = a(T — T C )M 2 + bM 4 We are only interested in the case of vanishing external 
magnetic field. Above T c , the lowest free energy state has zero net magnetization, while 
below T c the material acquires a magnetization M — ±[a|T — T c \/2b] 1 ^ 2 . 

It is also well known that mean-field theory breaks down close to T c ; for example, al- 
though mean-field theory correctly predicts that the correlation length diverges as £ ~ 
|(T — T c )/T c \~ u , it gives the value of the critical exponent v = 0.5, while numerical sim- 
ulations find v = 0.630(2) ||. In order to handle the infrared divergences that appear near 
T c , the renormalization group is used to relate a given theory to an equivalent theory with 
smaller correlation lengths and thus better behaved in the infrared. With the ^-expansion, 
one works in 4 — e dimensions and finds a fixed point of order e of the renormalization group 
equations, taking the limit e — > 1 in the end. To second order in e, v — | + + j^z 2 — 0.63, 
a remarkable result ||. 

In the present work, we would like to follow a somewhat orthogonal approach to study a 
closely related question. Consider again an Ising ferromagnet in the absence of an external 
magnetic field, which is well below its critical temperature and has been prepared with all 
spins initially pointing in one direction. In the thermodynamic limit, this is the broken- 
symmetric state. As the temperature is increased, thermal fluctuations will flip groups 
of spins, so that the absolute value of the magnetization will start decreasing. At some 
temperature T* < T c , \M\ — > 0. (This is another way of expressing the breakdown of mean- 
field theory.) The question we would like to address is the following: At a given temperature 
below T c , what is the fraction of the volume occupied by each of the two possible phases of 
the system as a function of the parameters of the potential? Can we provide a method for 
computing this fraction which somehow encompasses the breakdown of mean-field theory, 
without resorting to the renormalization group? 

We start by writing the homogeneous part of the free-energy density as 

y(0,r) = ^(r 2 -r 2 2 )0 2 -|T0 3 + ^ 4 . (i) 

It is easy to see that, with a simple field redefinition, we can rewrite V(4>, T) as a Ginzburg- 
Landau potential with a temperature-dependent magnetic field. At the critical temperature, 
this field vanishes and we recover the degenerate double-well potential. 

Introducing dimensionless variables x = a 1 ^ 2 T 2 x, t = a 1//2 T 2 t, X = aT^I^T^ (f>, and 
— T/T2, the Hamiltonian is, 



H[X] 1 



e e 



d 3 x 



2 1 
+ 2 



-*6X* + —X 4 
3 4 



1)X Z 

(2) 



where a = a 3//4 a, and A = a 1//2 A (henceforth we drop the tildes). At the critical tempera- 



ture 6 C = (l-2a 2 /9A)- 1 / 2 the two minima, at X = and X + = §f [l + yjl - 4A (1 - I/O 2 ) /a 2 
are degenerate. In what follows, we are only interested in the system at 8 C . 
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In a recent work, Borrill and Gleiser (BG) simulated the dynamics of the above system 
coupled to a Markovian thermal bath (This is why we wrote the potential as in Eq. |l[) 
In analogy with Ising ferromagnetism, it is useful to think of the system described by the 
Hamiltonian above as having two phases, the "0 phase", for X < X max , and the "+ phase", 
for X > X max , where X max is at the maximum of the double-well potential. The initial 
conditions were chosen so that the system started in the phase. Thus, V^/Vr{t = 0) = 1, 
where Vq is the total volume in the phase, Vr is the constant total volume, and V^/Vy(£) + 
V+/Vr{t) = 1. The coupling of this system to a thermal bath will induce the nucleation 
of domains of + phase. For a small free-energy barrier between the two phases, domains 
of the phase will also nucleate within domains of + phase, and so on, resulting in a 
very complicated domain structure. It is necessary to distinguish between two possible 
kind of domains, connected and disconnected. Connected domains percolate throughout the 
volume, and thus cannot be surrounded by a domain of the opposite phase. Disconnected 
domains have finite volume and are always surrounded by the other phase. In general, 

T/i _ T/c _|_ yd 

V 0(+) - V 0(+) + v o(+)- 

Keeping the system at 6 C and fixing a = 0.065 (this seemingly ad hoc choice was inspired 
by the electroweak effective potential, although any other value would do), BG measured the 
equilibrium fractions, fo = Vq/Vt and /+ = V+/Vt, as a function of the coupling A. They 
found that for A > A c ~ 0.025, fo = f+ = 0.5, that is, the two phases completely mix (see 
Fig. 1). In other words, for A > A c , the symmetry is restored, even though the mean-field 
potential still has a double-well shape. Clearly, for A > A c , mean-field theory breaks down, 
and the system is better described by an effective free-energy density with a single minimum 
at X max = 0. This situation is exactly analogous to an Ising ferromagnet for T >T*. 

In what follows, we will reproduce these results with a simple statistical model for the 
thermal fluctuations. Our approach is completely general, in that it can be easily adapted to 
other systems described by a similar double- well potential, i.e., for systems in the universality 
class of the Ising ferromagnet. We will assume that the large- amplitude fluctuations from the 
0(+) phase into the +(0) phase can be modelled by spherically-symmetric subcritical bubbles 
of Gaussian shape of a given radius and amplitude (For previous treatments of subcritical 
bubbles see Ref. ||.), 

Mr) = M- r2,R2 , Mr) = <Pc (l - e- r2 l R2 ) , (3) 

where R is the radial size of the configuration, and <ft c is the value of the field amplitude at 
the bubble's core, away from (and into) the phase. For these configurations to interpolate 
between the two phases in the system, 4> c > max . With this ansatz, the free energy of the 
fluctuations assumes the general form, F SC (R, <fi c , T) = &<^i?+C/3(0 c , T)R 3 , where b, and c are 
numerical constants and f3((p c ,T) is a polynomial which depends on the particular potential 
used. We will further assume that the nucleation rate for these configurations is obtained 
from a Gibbs distribution, G(0, R) = Aexp [— F sc /T], where A is a constant. The nucleation 
rate per unit volume is then r = / Gdifi dR. 

The number density of bubbles of, say, the + phase with radii between R and R + SR 
and amplitudes between <p and <fi + 5<p at time t + 5t is 

n + (R + 5R,(f) + 5<f>,t + 5t) -n + (R,4> + 5(f),t + 5t) 
- [n+(R + 6R,<l>,t + 6t)-n + (R,(f>,t + 6t)] . 
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Changes in the number density are generated by three different processes: i) bubbles can 
shrink into and out of this interval. [We assume the time dependence of the amplitude is 
closely related to the radial time-dependence, as recent numerical studies have demonstrated 
|J.]; ii) bubbles can be thermally nucleated into this interval; hi) bubbles can be thermally 
destroyed out of this interval by inverse nucleation. Expanding to first order in 5R, 5(f), and 
5t, we obtain a Boltzmann equation for the bubble distribution function f^(R,<fi,t), 

dt - ~^m + V T G ^ + 

V? 

TT G +^o , (4) 



Vi 



T 



where the bubble density distribution function for domains of the + phase, [hence the su- 
perscript d, for disconnected], is defined as 

and G r (+)-*+(o) is the Gibbs distribution for the nucleation rate per unit volume of bubbles 
of the +(0) phase within the 0(+) phase. Note that we have written the Boltzmann equation 
to be consistent with the initial conditions used in the simulations, so that only disconnected 
domains contribute to the fraction of the volume in the + phase. It is straightforward to 
adapt the equation to different initial conditions. 

In order to proceed, we note that the total fraction of volume in the + phase can be 
written as 

7 ~ / / — — —±d<f)dR . (6) 




3 J d<pdR 

For a degenerate double-well, G+^o = Go^+ = G, and the Boltzmann equation can be 
written as, 

. _ H M + (1 _ 27)G . 

Imposing the physical condition f+(R — > oo, 0, t) — > 0, we can solve for the equilibrium 
distribution function (/? = 0), and use it to compute the equilibrium fraction of the volume 
in the + phase, 7 eq , which is the quantity measured in the BG simulation. The general 
solution is, 



7eq 



where, 



-^(0max; Rtdivl) 
1 2/(0 m ax? -Rmin) 



J(0max, #min) = fT / / / dR'dRd(f) 



|^| JR -'-Rmin "tymax 

An 

x—R 3 exp[-F sc (R'A)} , (9) 
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and max is at the maximum of the double well potential, and the minimum radius is taken 
to be the lattice spacing in the BG simulation, i.e. R m i n = 1, which sets the coarse-graining 
scale. An analytical expression for 7(0 max , Rmm) can be obtained if we write F SC (R', <f>) ~ 
b(p 2 R'. This is a good approximation in the case that fluctuations are small enough that the 
volume term is sub-dominant. 

In Fig. 1, we show 7 eq as a function of A. The dots are the results of BG, while the 
curves are the results of the integration of the Boltzmann equation, for different values of 
the single parameter It is clear that for A < A c , we obtain an excellent fit to the data 

with A/\v\ = 60. 

For larger values of A, our method underestimates the fraction of volume in the + phase. 
This can be understood by noting that our kinetic description does not include possible terms 
which account for the coalescence of nearby domains. As fluctuations become more probable, 
these terms will play an increasingly important role. In fact, it is easy to understand how the 
onset of domain instability is intimately related to percolation of the + phase. (An excellent 
introduction to percolation theory can be found in Ref. 0.) 

Consider a large spherical domain of the + phase of radius R. There are three ways by 
which the volume V of this domain can change; i) it may shrink by surface tension with a 
velocity v ; ii) a small bubble of the + phase may nucleate just outside it; iii) a small bubble 
of the phase may nucleate inside it. Assuming that the bubble of the phase nucleates just 
inside the large domain (see Fig. 2), we can write an approximate equation for the change 
in the volume V, 

^ = -vAnR* + (T ^ + AV) ^r 3 - (T + ^ AV) ^r' 3 , (10) 
at 3 3 

where ro(+)->+(o)AV( ) is the nucleation rate for a bubble of the +(0)-phase of radius in 
the neighborhood of the domain wall. Assuming for simplicity that r = r', and recognizing 

that r_| >o(r) = r ^+(r) for a degenerate double-well potential, the condition for domain 

instability, ^ > 0, becomes, 

rr 4 >A^),. (ii) 

On the other hand, in order for the + phase to percolate, 7 eq > p c , where p c is the critical 
percolation probability. Using Boltzmann's equation, we can approximately write 

Q 47T A , . 

7eq = — — , a ^ — rr 4 . 12 
/q 1 + 20' y 3v K J 

Thus, for percolation, we obtain the inequality, 

rv > ^^— v . (13) 

For a simple cubic lattice, p c = 0.311. Note that this is remarkably close to the value of 7 eq 
for A c , the point where the kinetic description breaks down (see Fig. 1). 

Moreover, comparing the two inequalities of Eqs. [TT] and [13], it is clear that they are 
satisfied at similar values of Tr 4 . In particular, writing for simplicity p c = 1/3, they are 
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equal for R = 2r. This simple argument strongly suggests that the onset of percolation is 
dynamically driven by the nucleation of small bubbles in the neighborhood of large domains. 

As stated before, once the barrier between the phases is small enough, the domain wall 
becomes unstable and nucleation of small bubbles is the dominant process for changing a 
volume of space from one phase to another. The equations determining the fraction of 
volume in each phase can then be approximated as 



(yr 3 iw) V (t) - (f r 3 r + ^ ) V + (t) 



dV + fA-K \ _ /47T 

V Q (t) = V w -V+(t), (14) 



where r is the average radius of a nucleated bubble. The solution to these equations is 

V + (t)_l , fV + (0) _ 1\ e _ {87r r r 3 /3 ). (15) 



V tot 2 V V tot 

Notice that the equilibrium value V + (oo)/Vt t = 7eq = 1/2, independently of V+(0). 

Therefore, once the domain wall becomes unstable, the system will equilibrate to 7 = 1/2 
in a time scale ~ (Tr 3 ) -1 , which, at percolation, is of order the light crossing time r (see 
Eq. [13]). That is, once the system percolates, the stable, equilibrium, mean value of the field 
is at max , which is the value exactly between the two minima. 

As shown eq. (|15D , even if the system is perturbed away from 7 = 1/2, it will quickly 
relax back to this value. This indicates that the percolation point is the point of symmetry 
restoration: the symmetry of the true (coarse-grained) effective potential has been restored, 
even though the mean field potential still describes a double-well potential. It is important 
to note, however, that even if the average value of the field is maX ; there are widespread 
large amplitude fluctuations with average volume set by the correlation length; the system 
is far from being locally homogeneous. 

We conclude that our method based on subcritical bubbles gives a quantitatively accurate 
description of the behavior of thermal fluctuations for models in the Ising universality class. 
It also provides a dynamical picture of symmetry restoration, and the breakdown of mean- 
field theory. Furthermore, the method relates the breakdown of mean-field theory to a critical 
value of a given parameter, which is easily calculable. For values of the parameter larger than 
the critical value, large domains become unstable to growth due to the nucleation of nearby 
bubbles, percolation ensues, and symmetry or, more accurately, complete phase mixing, is 
restored. 
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Figure 1. The fraction of the volume in the + phase. The dots are from the numerical 
simulations of BG, while the lines are the solutions of the Boltzmann equation for different 
values of the parameter 

Figure 2. Schematic of domain instability. The surface tension on bubble wall will tend 
to shrink the large bubble with wall velocity v. However, more bubbles of + phase will 
nucleate just outside the large bubble wall than bubbles of phase just inside the wall 
because AV > AV (nucleation rates of and + bubbles are equal). This will tend to make 
the large bubble grow. When the barrier between the phases is small enough, nucleation 
dominates over bubble shrinking, causing the wall to become unstable to rapid growth. 
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